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ABSTRACT 

Aims. This study aims to characterize linear polarization structures in LOFAR observations of the interstellar medium (ISM) in the 
3C196 field, one of the primary fields of the LOFAR-Epoch of Reionization key science project. 

Methods. We have used the high band antennas (HBA) of LOFAR to image this region and Rotation Measure (RM) synthesis to 
unravel the distribution of polarized structures in Faraday depth. 

Results. The brightness temperature of the detected Galactic emission is 5 - 15 K in polarized intensity and covers the range from 
-3 to +8 rad m“^ in Faraday depth. The most interesting morphological feature is a strikingly straight filament at a Faraday depth of 
+0.5 rad rcT^ running from north to south, right through the centre of the field and parallel to the Galactic plane. There is also an inter¬ 
esting system of linear depolarization canals conspicuous in an image showing the peaks of Faraday spectra. We used the Westerbork 
Synthesis Radio Telescope (WSRT) at 350 MFlz to image the same region. For the first time, we see some common morphology in 
the RM cubes made at 150 and 350 MHz. There is no indication of diffuse emission in total intensity in the interferometric data, in 
line with results at higher frequencies and previous LOFAR observations. Based on our results, we determined physical parameters 
of the ISM and proposed a simple model that may explain the observed distribution of the intervening magneto-ionic medium. 
Conclusions. The mean line-of-sight magnetic field component, 6||, is determined to be 0.3 ±0.1 pG and its spatial variation across 
the 3C196 field is 0.1 pG. The filamentary structure is probably an ionized filament in the ISM, located somewhere within the Local 
Bubble. This filamentary structure shows an excess in thermal electron density («eS|| > 6.2 cm^^pG) compared to its surroundings. 

Key words. ISM: general, magnetic fields, structure - radio continuum: ISM - techniques: interferometric, polarimetric 


1. Introduction 


The Galactic interstellar medium (ISM) is filled with cold, warm. 


and hot thermal gas (the multi-phase medium; e.g. McKee & Os- 
triker|19T7) in a mixture of ionized, neutral, atomic, and molec¬ 
ular components. In addition, the ISM is permeated with non- 
thermal plasma (mostly relativistic protons and electrons) per¬ 
vaded by a large-scale magnetic field. In general, the energy den¬ 
sity in the various components is approximately equal, but large 
deviations in energy balance occur in various situations. Obser¬ 
vations of diffuse polarized emission are valuable in exploring 
both the distribution and the properties of the ISM. Magnetic 
fields, relativistic electrons, and thermal electrons can be con¬ 
strained through the action of Faraday rotation, which changes 
the morphology of the observed emission. 

Studies of Galactic polarized emission are mostly done at 
high radio frequencies (> 1 GHz; for an overview see [Reich 
|2006| and references therein). At lower radio frequencies, these 
studies were performed mainly at 325 MHz with the Westerbork 
Synthesis Radio Telescope (WSRT). These studies revealed a 
large number of unusually-shaped polarized, small-scale struc¬ 
tures, which have no counterpart in the total intensity (e.g. 


Wieringa et al.|1993 Haverkorn et al.|2003a|b| Schnitzeler et al. 


2009|l. They are attributed to a combination of Faraday rotation 
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and depolarization effects due to variable Faraday depth in the 
ISM. In many areas, linear structures dominate the morphology 
and have so far defied explanation. This is especially the case 
at higher Galactic latitudes where line-of-sight superposition ef¬ 
fects would not inherently randomize anisotropic structures. 

Faraday rotation is proportional to the square of the wave¬ 
length. Therefore, low radio frequency observations are very 
sensitive to small column densities of magnetized plasma in 
the ISM that are difficult to detect at higher radio frequencies. 
The wide frequency coverage and good angular resolution of the 
Low-Frequency Array (LOFAR; |van Haarlem et al.|2013| l, make 
LOFAR an excellent instrument for studying Galactic polarized 
emission Pacobelli et al.||^13[ jJelic et ^|2014] l. In combina¬ 
tion with the rotation measure (RM) synthesis technique ( jBren-j 
tjens & de Bruyn|2005| l, LOFAR observations allow us to study 
the relative distribution of synchrotron-emitting and Faraday- 
rotating regions at an exquisite resolution of ~ 1 rad m“^ in 
Faraday depth. In addition, the high angular resolution available 
in LOFAR allows us to study the medium at a resolution hardly 
affected by beam depolarization. 

Recent observations with LOFAR revealed diffuse polariza¬ 
tion in several fields at high Galactic latitudes with surprisingly 
high brightness temperature. For example, polarized emission in 
the ELAIS-Nl field is ~ 4 K at 150 MHz Pelic et al.||20T^ . 
This is much greater than was anticipated on the basis of earlier 
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Table 1. An overview of observational parameters of the five LOFAR- 
HBA observations used in this study. We are using in our analysis only 
6h, symmetric around transit, of each 8h observation. 


Observation ID 

Start Time [UTC] 

L79324 

06-Dec-2012 22:41:05 

L80273 

12-Dec-2012 22:17:30 

L80508 

16-Dec-2012 22:01:46 

L80897 

21-Dec-2012 22:42:46 

L192832 

15-Dec-2013 23:06:40 

Phase centre (J2000.0) 

0843”36.07% h-48°13'02.58" 

Frequency range 

115-189 MHz 

Spectral resolution 

3.2 kHz 

Integration time 

2 s 

Observing time 

8(6) hours in 2012 (2013) 


WSRT (e.g. IBemardi et al]|2009l |2m0| |Pizzo et al.|[20TT]l and 
Giant Metrewave Radio Telescope (GMRT; e.g. Pen et al.|2009| l 
observations in the same frequency band. Early results from the 
Murchison Widefield Array (MWA) also indicated much higher 
levels of observed polarization ( |Bernardi et al. 12013) 1. In this pa¬ 
per, we present LOFAR observations of diffuse polarized emis¬ 
sion detected in the 3C196 held, one of the primary helds of the 
LOFAR-Epoch of Reionization key science project. 

The 3C196 held is centred on the very bright quasar 
3C196, in a cold region towards the Galactic anticentre (1=171°, 
b=H-33°). It has been observed previously with the Low Fre¬ 
quency Front Ends (LFFE; 138 - 156 MHz) on the WSRT radio 
telescope ( |Bernardi et al.|20T0) l. These data revealed only faint 
patchy polarized emission, restricted to Faraday depths smaller 
than 4 rad m“^. The surface brightness of this emission was close 
to the thermal noise on angular scales smaller than 10'. On scales 
greater than 30', emission has an rms value of 0.7 K. The reso¬ 
lution in Faraday depth space of the WSRT data was three times 
poorer than that of the LOFAR observations presented in this 
paper, which have three times wider frequency coverage. 

The paper is organized as follows. In Section 2 we give an 
overview of the observational set-up and the data reduction. The 
widefield images of the 3C196 field in total intensity and polar¬ 
ization are presented in Section 3. In Section 4 we discuss the 
hdelity of the data. In Section 5 we discuss properties of the de¬ 
tected diffuse polarized structures and give a possible theoretical 
interpretation, while analysis of the observed discrete polarized 
sources will be presented in a separate paper. We summarize and 
conclude in Section 6. 


2. Observations and data reduction 

The 3C196 field was observed multiple times during winter 
2012/13 and 2013/14 with the LOFAR High Band Antennas 
(HBA). During this period we accumulated ~ 450 hours of data. 
For this paper we have fully processed hve syntheses with dura¬ 
tions of 8h (2012) and 6h (2013), as a representative sample of 
the observations. An overview of these five observations is given 
in Table [1] 

The array was used in the so-called HBA DUAL INNER 
conhguration ( |van Haarlem et ar]|2013| l, consisting of 48 core 
stations (CS) and 14 remote stations (RS), which were tapered 
to have the same size and shape as the CS. Data were taken in 
the frequency range from 115 MHz to 189 MHz, divided into 
380 sub-bands of width 195.3125 kHz. Each sub-band was fur¬ 
ther divided in 64 channels. The correlator integration time was 
2 s. All four correlation products between pairs of orthogonal 
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Fig. 1. Frequency-averaged Stokes I image of the 3C196 region, 
based on the four nights of data taken in December 2012. The im¬ 
age is 8.3° X 8.3° in size, with a PSF of 45", and the noise level 
is 170 pJyPSF^*. The nine bright sources towards which we per¬ 
formed direction-dependent calibration are indicated with red circles. 
The 3C196 source is subtracted from direction-independent calibrated 
visibilities, while eight other sources are not. A location of the discov¬ 
ered pulsar J081558-I-461155 (see Sec. |3.3[ > is shown with a blue as¬ 
terisk. The image is de-convolved using the CLEAN algorithm imple¬ 
mented in CASA. 
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dipoles were recorded. The total time per observation was 8 h in 
2012 and 6 h in 2013. All observations were symmetric around 
transit and were taken during the nighttime. The uv coverage was 
fully sampled up to baselines of about 800 wavelengths. To com¬ 
pare observations taken in different observing cycles, we used 6h 
symmetric around transit of each 8h observation. 

The processing of the data was done on a CPU/GPljQcluster 
dedicated to the LOFAR-EoR project, located at the University 
of Groningen, the Netherlands. During the processing we used 
the LOFAR-EoR Diagnostic Database (LEDDB), which is used 
for the storage and management of the LOFAR-EoR observa¬ 
tions ( |Martinez-Rubi et al.|2013| l. 


2.1. Initial processing and calibration 


The first step in our initial processing is the automatic flagging of 
radio frequency interference (RFI) using AOFlagger ( jOffringa] 
et al.pOlO 2012| l. In the frequency range from 115-177 MHz 
on average only ~ 3% of our data is flagged. This percentage, 
however, is much higher (>40%) for frequencies above 177 MHz 
where signals from Digital Audio Broadcasting (DAB) corrupt 
the data in a few 1.75 MHz wide bands. For an overview of the 
LOFAR RFI-environment we refer to Offringa et al. (2013|l. 


After flagging, the data are averaged to reduce the data vol¬ 
ume for further processing. The bandwidth and time smearing 
resulting from this averaging only affects the longer baselines, 
which are not used in this study. To remove edge effects from the 
polyphase hlter the 4 edge channels from the initial 64 channels 
are excluded in the averaging process. The resulting data set has 
a spectral resolution of 183 kHz (henceforth called a sub-band) 
and a time resolution of 10 s. 


* CPU: Central Processing Unit; GPU: Graphics Processing Unit 
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V. Jelic et al.: Linear polarization structures in the 3C196 field 



Time [h] 

Fig. 2. The estimated RM variations in the ionosphere using the Global 
Ionospheric Maps of CODE. The RM are given as a function of observ¬ 
ing time and for each night. 


tracted from the residual visibilities with the directional gains 
derived for their directions. 


2.2. Faraday rotation in the Earth’s ionosphere 

Ionospheric Faraday rotation is a time- and direction-dependent 
propagation effect, whose size is proportional to the total elec¬ 
tron content (TEC) of plasma in the ionosphere. If variability 
in Faraday rotation happens on a timescale that is smaller than 
the total integration time of an observation, the observed polar¬ 
ized emission is partially decorrelated. To correct for this effect 
we predict the RM variations in the ionosphere using the Global 
Ionospheric Maps of the Centre for Orbit Determination in Eu¬ 
rope (CODiQ and the World Magnetic Model (WMi yp]). A full 
description (rf this procedure is given in|Jelic et al.|(|2014| sec. 
3.2). 

Figure shows the estimated RM variations in the iono¬ 
sphere that we apply to the data. The curves are given for dif¬ 
ferent nights as a function of observing time. The RM variations 
are the largest during the L79324 observation. This observation 
also shows the largest absolute RM values. We emphasize that 
at 150 MHz, a RM variation of 0.2 rad m“^ rotates the plane of 
polarization by 45°. 


The calibration proceeds in two steps. First, we perform 
a direction-independent calibration with the Black Board 
Selfcal (BBS) package (Pandey et al. 2009 | l. Each sub-band 
is calibrated separately. The source 3C196, which is located in 
the pointing and phase centre, dominates the visibilities on all 
baselines, making this step of the calibration straightforward. 
We use a four-component model for 3C196, which accurately 
represents the high-resolution structure of the source (Pandey et 
al., in prep.). The total flux of 3C196 (84 ± 2 Jy at 150 MHz) 
is set by the broadband spectral model given by |Scaife & Heald] 
( 2012| l. The BBS step removes the clock and ionospheric phase 
errors, and sets the frequency-dependent intensity as well as the 
astrometric reference frame for the field. The BBS solves for all 
four complex elements of the station Jones matrices taking the 
changing location of 3C196 within the dipole element beam into 
account. Because the station phased array centre tracks 3C196, 
the varying station gain is automatically accounted for in this 
process. When applying the gains we also correct for the varying 
parallactic angle. This leads to a very stable and small polarimet- 
ric response over the inner part of station beam. 

The second step is more involved. There are at least eight 
additional sources with integrated apparent flux densities in ex¬ 
cess of 2 Jy (at 150 MHz) within 2.5° of 3C196. To remove 
their time-variable polarized response, together with one of the 
3C196, requires direction-dependent calibration. This is done 
using S AGECal ([Yatawatta et al. 2008 [ Kazemi et al. 2011 
|2013a|bHKazemi & Yatawatta|2013 i in nine directions. This cor¬ 
rects for time- and frequency-dependent errors due to the station 
beam and the slight position-dependent ionospheric delays. Each 
direction is associated with one source. The solution interval is 
20 minutes for the eight bright sources and 2 minutes for 3C196. 
This interval is sufficient to remove tfiese sources down to tfie 
confusion noise on tfie short core baselines. During SAGECal we 
are only using baselines > 800 wavelengths to minimize sup¬ 
pression of the ubiquitous extended diffuse (polarized) emission, 
which is not included in the sky model and mostly appears at 
baselines < 800 wavelengths. All nine sources are unpolarized 
down to a fraction of a percent. These nine sources (3C196, and 
eight bright, discrete sources; see Fig.[T]| were subsequently sub- 


2 . 3 . Imaging and RM synthesis 


If not stated otherwise, images we present are produced using 
the excon imager ( | Yatawatta||20 1 4[ http;//exconimager.sf.net/). 
To analyse the polarization of diffuse emission, which mostly 
appears on spatial scales greater than a few arcmin, we make 
lower resolution images in all Stokes parameters (IQUV). These 
images are produced using baselines between 10 and 800 wave¬ 
lengths, providing a frequency independent resolution of about 
3 arcmin. We use robust (Briggs) weighting with robustness pa¬ 
rameter equal to 0. 


To study the linearly polarized emission as a function of 
Faraday depth (O), we apply RM synthesis (Brentjens & de 


|Bruyn|2005| l to the Stokes Q,U images of 310 sub-bands, which 
have comparable noise levels. These sub-bands cover the fre¬ 
quencies between 115 MHz and 175 MHz. We first syntfie- 
size a cube over a wide range in Faraday depth to determine 
where polarized emission could be detected. The final cube cov¬ 
ers a Faraday depth range from -25 rad m“^ to -1-25 rad m“^ in 
0.25 rad m"^ steps. The resolution in Faraday depth space, de¬ 
fined by a width of a RMSF (rotation measure spread function), 
is 50 = 0.9 rad uT^, while the largest Faraday structure that 
can be resolved is AO = 1.1 rad m“^. Since the resolution is 
comparable to the maximum detectable scale, we can only de- 


tect Faraday thin structures (4^ AO 
[2005] l. 


Brentjens & de Bruyn 


To determine the brightness temperature of emission and its 
angular power spectrum, we correct the RM cube for the relative 
variation of the element beam pattern across the field of view 
and for the primary beam. These corrections are applied using 
the AWimager ( |Tasse et al.|2013| l. We have not attempted to de¬ 
convolve the RM cube for the effects of the side lobes of the 
RMSF. The sidelobes are small and the S/N is generally so low 
that this would have had very little effect on the images. 


^ http://aiuws.unibe.ch/ionosphere 
^ http://www.ngdc.noaa.gov/geomag/WMM 
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L80508 BBS, PI: $=0.0 rad/m^ L80508 SAGECal, PI: $ = 0.0 rad/m^ 



128 126 124 122 120 118 128 126 124 122 120 118 

a (J2000) a (J2000) 


Fig. 3. Two images of the 3C196 region in polarized intensity at Faraday depth 0 rad m ^ after performing direction-independent calibration 
(image on the left) and direction-dependent calibration and source subtraction (image on the right). Images are 8.3° x 8.3° in size, with a PSF 
of 3.6' X 3.9'. The 3C196 source is also subtracted from direction-independent calibrated visibilities. Eight bright sources, used for direction- 
dependent calibration and source subtraction, are indicated with red circles. The images have not been corrected for polarization noise bias nor for 
Faraday rotation in the Earth’s ionosphere. 


3. Observational results 

3.1. Direction-dependent caiibration: instrumentaiiy poiarized 
sources 


The beam pattern of a LOFAR-HBA station can be described as 
the product of an element beam pattern and the array beam pat¬ 
tern of the station ( |van Haarlem et al.|2013| l. The element beam 
pattern is strongly polarized. Its polarization response is related 
to the projection of two beam patterns of two orthogonal dipoles 
on the sky and the changing parallactic angle. As a result of this, 
spurious polarization is produced. We correct the data for the 
beam pattern during calibration to first order, using BBS. The 
data are corrected for both the array and the element beam gain 
at the centre of the image. However, the relative variation of the 
element beam pattern across the field of view, as well as the tem¬ 
poral changes, are still be present in the data (for more details, 
see [Asad et al. [20T5] l. The direction-dependent calibration and 
source subtraction can correct for these effects to some extent. 

The results of direction-dependent calibration and source 
subtraction are illustrated in Fig. We show images in polar¬ 
ized intensity before and after performing direction-dependent 
calibration and source subtraction. Images are shown at a Fara¬ 
day depth of 0 rad m“^. The instrumental polarization leakage, 
after calibrating on 3C196, is to first order independent of fre¬ 
quency, and therefore all instrumentaiiy polarized signals are lo¬ 
cated around a Faraday depth of 0 rad m“^. 

All sources included in the sky model are removed success¬ 
fully. Their side lobes are significantly suppressed and they are 
not visible in the images. The morphology of diffuse emission 
is preserved. A flux scale of diffuse emission agrees on average 
over all directions within 5 - 10% with the data calibrated using 
only BBS (direction-independent calibration). 


3.2. Diffuse poiarized emission 

A series of widefield images of the 3C196 field in both polar¬ 
ized intensity and Stokes Q,U are presented in Fig. The im- 
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ages are 8.3° x 8.3° in size with a PSF (point spread function) of 
3.6' X 3.9'. The noise level is 71 piy PSF'' RMSF'* in polarized 
intensity and 108 pJy PSF * RMSF * in Stokes Q,U. The im¬ 
ages are given at Faraday depths of-1.5, 0, -i-l, -1-1.5, h- 2, -1-2.5, 
-1-3.0, and -1-3.5 rad m'^ to emphasize the various detected struc¬ 
tures in linear polarization. 

Diffuse emission (< 3 mJy PSF'* RMSF'*) is detected over 
a range of Faraday depths from -3 to 4-8 rad m'^. A triangu¬ 
lar feature in the south-west side of the primary beam (marked 
with C in Fig.|^ dominates at negative Faraday depths. At pos¬ 
itive Faraday depths, emission first appears as two extended 
north-south structures (denoted with A and B in Fig. |^. A 
striking filamentary structure A, running right across 3C196, 
is at least 4° long and 7.5' wide. Its mean surface bright¬ 
ness is 3.2 mJy PSF'* RMSF'* and it peaks at Faraday depth 
~ 0.5 rad m'^. The two-legged puppet structure B has the mean 
surface brightness of 2.5 mJy PSF * RMSF'*, peaking at Fara¬ 
day depth ~ 0.25 rad m'^. Towards the higher Faraday depths, 
emission fills up the east side of the primary beam. The mean 
surface brightness of this emission is 1.6 mJy PSF * RMSF'*. 


A few representative Faraday spectra of detected features are 
shown in Fig. |5 The largest Faraday structure that can be re¬ 
solved is 1.1 rad m'^ wide. Information contained in the Fara¬ 
day spectra is summarized in Fig. ^ We follow Schnitzeler et al. 
P009| l and show a map of the highest peak of the Faraday depth 
spectrum at each spatial pixel and a map of the RM value of each 
peak. The highest peaks in the Faraday spectra have on average a 
surface brightness of ~ 2.5 mJy PSF'* RMSF'*, located mostly 
at positive Faraday depths. 


While the polarized intensity shows only the amplitude of 
polarized emission, images in Stokes Q,U also reflect the mor¬ 
phology of polarization angle. There are a number of striking 
patterns and features in Stokes Q,U (see Fig. |^. Most of them 
are outlined by regions that show emission in Stokes Q,U of op¬ 
posite sign. They have a typical width of 5' - 10'. 
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L80508, PI: $ = + 1.5 rod/m^ 


L80508, Q: $= + 1.5 rod/m^ 


L80508, U: $ = +1.5 rad/m^ 
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Fig. 4. Widefield images of the 3C196 field in polarized intensity (PI), Stokes Q, and Stokes U at Faraday depths of -1.5, -0.0, -l-l.O, -1-1.5, -1-2.0, 
-1-2.5, -^3.0, and -t3.5 rad nr^. Images are 8.3° x 8.3° in size with a PSF of 3.6' x 3.9'. The noise level is 71 yujy PSF“* RMSF^' in polarized 
intensity and 108 pJy PSF^' RMSF^' in Stokes Q,U. The images in polarized intensity have not been corrected for polarization noise bias nor for 
the primary beam. 
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L80508, Q: $ = + 2.0 rod/m^ 



128 126 124 122 120 118 

a {J2000) 

L80508, Q: $ = + 2.5 rod/m^ 


L80508, U: $ = + 2.0 rad/m^ 
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L80508, U: $ = + 2.5 rad/m^ 



L80508, PI; $ = + 3.0 rad/m^ 


L80508, Q: $ = + 3.0 rad/m^ 


L80508, U; $ = + 3.0 rad/m^ 



L80508, PI: $ = + 6.5 rod 


i/m^ 
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L80508, Q: $= + 6.5 rod/ 
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Fig. 4. Continued. 


3.3. Polarized sources 

In addition to the widespread diffuse emission, we identified 
~ 15 discrete polarized sources in the 3C196 field. All but one 


of these sources are (giant) radio galaxies with typical polariza¬ 
tion fractions of a few percent, and a detailed analysis of their 
properties will be presented in de Bruyn et al. (in prep). One 
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V. Jelic et al.: Linear polarization structures in the 3C196 field 


L80508 




Fig. 5. Faraday spectra of a few representative lines-of-sight through the 
RM cube of polarized intensity. We have not attempted to deconvolve 
the spectra for the effects of the side lobes of the RMSF. 


source, J081558+461155 (see Fig. [T]&|^, has a steep spectrum, 
is spatially unresolved and shows about 50% polarization. These 
properties are reminiscent of those of pulsars. 

A follow-up observation of this field, in Director’s Dis¬ 
cretionary Time, with the LOFAR Full HBA Core in tied- 
array mode ( [van Haarlem et ar]|2013| l indeed revealed a pul¬ 
sar with a period of 434 ms and a dispersion measure (DM) of 
11.28 pc cm“^ (private communication, J. W. T. Hessels and V. 
Kondratiev; see the LOFAR Tied-Array All-Sky Survey Discov¬ 
eries. 


The RM of the pulsar is h- 2.7 ±0.1 rad m“^, as determined 
from our RM cubes. In the direction of the pulsar the NE2001 
model of the Galactic distribution of thermal electrons ( |Cordes| 
|& Lazio] |2002| |2003 1 predicts a thermal electron density of 
(we) = 0.028 cm“^. Given its DM, we estimate its distance to 
be ~ 400 pc. Given the sparsity of sources in this direction, how¬ 
ever, both the density and the distance derived from the model 
must be considered uncertain. 


4. Fidelity of the data 

In this section, we make qualitative and quantitative comparisons 
between five LOFAR observations. The analysis is performed on 
RM cubes in polarized intensity, Stokes Q and U. Table gives 
an overview of the noise in the RM cubes for different nights. 
Variation of the noise from night to night is small. The L79324 
observation shows the largest variations in RM and has the high¬ 
est noise. L80273 has the lowest noise. The ratio of the rms val¬ 
ues between the polarized intensity and Stokes Q,U is ~ 0.65. 
This implies consistent normally distributed noise in Stokes Q,U. 

We then compare maps of the highest peak of the Faraday 
depth spectrum at each spatial pixel for different nights. An ex¬ 
ample of this kind of map is given in Fig. For each pair of 
maps, we calculate the Pearson correlation coefficient using only 
the inner 5° x 5° of the images. We find that the peak fluxes 

http://www.astron.nl/lotaas/ 


Table 2. The noise in the RM cubes of different observations. The val¬ 
ues are given for polarized intensity (PI) and Stokes Q,U. 


Observation 

PI 

Q 

U 

ID 

[juJy PSP * RMSE-i] 

L79324 

83 

123 

127 

L80273 

67 

101 

102 

L80508 

70 

111 

105 

L80897 

75 

118 

114 

L192832 

73 

111 

115 


agree well, with a correlation of 96 + 1% from night to night. 
The L80273 and L80508 observations show the highest correla¬ 
tion, of 98%. This shows that LOFAR performance from night 
to night is very stable and that we are able to correct for Faraday 
rotation in the ionosphere to first order. The small difference in 
the peak fluxes between the observations is mainly caused by the 
low time resolution of the CODE data (2 h). 

To check the alignment of the RM cubes in Earaday depth, 
we calculate the cross-correlation between two observations as 
a function of the lag in Earaday depth. In this comparison 
we take only spectra that have a peak flux in polarized inten¬ 
sity > 1.5 mJy PSE^' RMSE^' (at least two times above the 
noise). We then average the cross-correlation functions to de¬ 
termine their common peak by fitting a parabola. The posi¬ 
tion of this peak gives the misalignment of the two RM cubes. 
On average we find that the five RM cubes are misaligned by 
0.10 ± 0.08 rad uT^. The best alignment is between L79324 and 
L80508 (0.02 ± 0.01 rad m“^) and the worst is between L80508 
and L192832 (0.14±0.01 rad m“^). A misalignment between dif¬ 
ferent observations can be attributed to systematic errors in the 
CODE data, resulting in a bias in TEC maps from day to day. 
This bias leads to uncertenties of ~ 1 TEC unit, corresponding 
to an error in RM around 0.1 rad m“^. 

To get a better signal-to-noise ratio in the RM cubes, one 
can average the cubes from different nights. However, the RM 
cubes should be aligned in Earaday depth to high accuracy. Even 
a small misalignment decorrelates the signal to some extent. The 
detected diffuse emission is ~ 25 x above the noise, so an ad¬ 
ditional improvement of Vs does not make an impact on our 
results and our current conclusions. 


Einally, we compare our LOEAR observations with obser¬ 
vations of the same field by Bernard! et al. (2010 1 in the same 
frequency range, obtained using the LEEE on the WSRT radio 
telescope. Emission detected with the WSRT-LEEE is fainter and 
its morphology appears much simpler (see Eigs. 9, 10, & 11 in 
[Bernard! et al.|20f0| l than the morphology of emission detected 
with LOEAR. There is just faint patchy emission, without any 
more distinctive morphological features. 


The difference between the WSRT-LEEE and LOEAR ob¬ 
servations can be attributed to the three times better resolution 
of LOEAR in Earaday depth. Multiple Earaday-thin structures 
(Affi < 50wsrt-lffe) detected in the LOEAR observations decor¬ 
relate when we observe them with a broader RMSE if their emis¬ 
sion shows different polarization angle. Moreover, WSRT-LEEE 
observations were showing some instrumental artefacts that con¬ 
taminated RM cubes ( jBemardi et al.|20T0| l. Given these restric¬ 
tions a more quantitative comparison between the LOEAR and 
WSRT-LEEE observations is not very meaningful and we restrict 
our analysis to the superior quality LOEAR data. 
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Fig. 6. The upper image shows the highest peak of the Faraday depth spectrum in polarization at each spatial pixel, while a map of the RM value 
of each peak, shown in red-blue colour scale, is on the bottom. Yellow contours are showing Ha intensity on a one degree scale ( [FinkbeinCT 2003 1 
[Haffner et al.|200^ , given in units of Rayleighs (see Sec. [5)1. A white box indicates the boundaries of the WSRT 350MHz image in Fig. 10| A 
location of the discovered pulsar J081558-I-461155 (see Sec.|3.3[> is shown with a yellow asterisk. 
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Fig. 7. Angular power spectrum of integrated polarized emission. The 
best power-law fit (solid line) to the power spectrum (“C/ - Z”) between 
I = 150 and Z = 2700 gives a slope of = -1.65 ± 0.03. Dashed lines 
show a range of data points used for the power-law fit. 


5. Discussion 


5.1. General properties of detected diffuse polarized 

emission and the underlying magnetic field component 


The diffuse emission presented in Sec. |3.2| is of Galactic origin. 
To analyse its properties, we follow our previous work (Jelic 
|et al.|[20T^ and compute the total polarized intensity at each 
pixel by integrating the polarized intensity RM cube along Fara¬ 
day depth. The brightness temperature of the resulting integrated 
emission varies from 5 to 15 K. This is comparable to the bright¬ 
ness temperature of emission detected in the ELAIS-Nl field 
( Jelic et ar]|2014| l, also a field at high Galactic latitude. Using 
only the inner 3° x 3° of the field, we estimate the angular power 
spectrum of the integrated polarized emission, C;. In the calcu¬ 
lation we assume the “flat-sky” approximation (e.g. White et al. 


1999 


trans 


, i.e. k^Pk - for Z = 2nk, where Pk is the 2D Fourier 


brm of an image and Z = 18O°/0[°] is a spherical harmonic 
multipole. Figure]^ shows the resulting averaged dimensionless 
angular power spectrum. We then fit a power law to “C/ - /” be¬ 
tween Z = 150 and Z = 2700. The best fit (solid line in Fig. 0 
gives a slope of Pi - -1.64 + 0.05. 


In total intensity a slope of the synchrotron emission is usu- 

ally steeper (-3.0 < Pi < -2.5 e.g. 

Baccigalupi et al. 2001 Tucci 

et al.||2002l |Bernardi et al.||2009| 

Ghosh et al.||20121l. A differ- 


ence between the two slopes can be explained by redistribution 
of power between different angular scales of polarized emission, 
mostly from large scales to small scales, as a result of small- 
scale structure of intervening Faraday-rotating medium. At LO- 
FAR frequencies even a very small RM of the medium leaves a 
signature in the power spectrum. 

On the same angular scales, 100 < Z < 3000, and in the 
same frequency range as our observations, but in an area located 
in the Galactic plane (Fan region), Bemardi et al.|([2009 1 found 


a slope for polarized emission of Pi - -1.65 + 0.15, which is 
very similar to the 3C196 field. The morphology of emission in 
the two fields is different, yet integrated emission along Faraday 
depth of both fields shows a very similar power spectrum slope. 


This indicates that in both fields redistribution of power through 
Faraday rotation is similar. 

In the current data we have not detected diffuse emission in 
total intensity. This can be attributed to a lack of short baselines 
in our observations. We are only sensitive to structures on scales 
smaller than two degrees, meaning that synchrotron emission in 
total intensity appears mostly on scales larger than two degrees. 
This is in line with the discussed steeper slope of the observed 
synchrotron emission in total intensity, and with previous obser¬ 
vations with LOFAR Pelic et al.|2014| l and WSRT (e.g. [Bernardi| 
|et al.|[20T0| ) at high Galactic latitudes. Unfortunately, a lack of 
detected emission in total intensity prevents us from making a 
more detailed study of Galactic interstellar turbulence through 
fluctuations in synchrotron emission, as was done by |Iacobelli| 
|et al.| ( |20T3] l for a low Galactic latitude field. 

In near future the AARTFAACj^(Amsterdam - ASTRON Ra¬ 
dio Transients Facility and Analysis Centre) project might help 
to mitigate a lack of a very short baselines in our current LOFAR 
observations. Six stations at the heart of LOFAR will be used as 
a 288-element array. This will then provide almost full uv cov¬ 
erage at short baselines and will enable us to map large-scale 
Galactic emission that is currently not accessible. 

To estimate the fraction of polarization, we divide the inte¬ 
grated polarized intensity by the 408 MHz total intensity map 
( Haslam et al.|198T 1982| l scaled to 160 MHz. The spectral in¬ 
dex between 45 and 408 MHz of Galactic synchrotron emission 
in this region is y6 = -2.5 ( Guzman et al.|20rT] l. If we scale the 
brightness temperature of ~ 23 K from 408 MHz to 160 MHz, 
we deduce a brightness temperature of 236 K in total emission. 
The observed polarization levels therefore imply a polarization 
of ^4%. 

A high degree of polarization indicates a regular large- 
scale magnetic field topology, coming from either uniform or 
anisotropically random magnetic fields. The uniform magnetic 
field component, following the spiral arms of our Galaxy, is al¬ 
most perpendicular to the line of sight at the direction of 3C196 
(for a description of Galactic magnetic field models, see Sun| 
|et al.||2008| and references therein). Its contribution to Fara¬ 
day rotation towards the 3C196 field is therefore expected to be 
small. 

This assumption is supported by the estimated mean line- 
of-sight magnetic field component in the direction of pulsar 
J081558H-461155 (see Sec. |3.3| l. Given the known relation be¬ 
tween the RM and DM of a pulsar (e.g. |Manchester|[l972|l, i.e. 


<B|I> _ RM[radm-2] 

“ 0.812 DM [pc cm-3]’ 


( 1 ) 


we get {Z?||) = 0.3 + 0.1 pG. A typical strength of magnetic field 
in the Galactic halo (thick disk) is ~ 3 pG, with an expected 
random component of ~ 1 pG ( Mao et al.|2010[ Haverkorn 2015| 
and references therein). 

Variations of the mean line-of-sight magnetic field compo¬ 
nent across the 3C196 field, can be estimated from vari¬ 

ations in the RM values of detected diffuse structures and in 
thermal electrons. We obtain the RM variations directly from 
the image in Fig. ctrm 4 rad m“^. A good tracer of warm 
ionized gas is the Hor line. From the composite all-sky Ha map 
|FinkbeinCT 2003| vl.l|^ we find the Ha intensity across the 


^ http://www.aartfaac.org 

^ This map is a composite of three separate surveys. The 3C196 field is 
covered only by the Wisconsin H-Alpha Mapper (WHAM;|Haffner et al.| 
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3C196 field to be 0.6 + 0.1 R (Rayleighs; 1 R = 10® (4;r)-^ pho¬ 


tons cm s ‘ sr 0- However, the given Ha intensity reflects the 
integrated emission over different velocities (|t;| < 100 km s“') 
and over the line-of-sight through the whole Galaxy, while the 
polarized structures detected with LOFAR are probably limited 
to distances < 1-1.5 kpc from us (see Sec. ig. Electron den¬ 
sities that we derive from Ha intensity should therefore be taken 
more as an upper limit. We do not know whether we see polar¬ 
ized emission all the way to the edge of the Galaxy. 

We convert the observed Ha intensity into emission measure 
(EM), according to ( [Reynolds |1 99 l| l. 


EM 


[pc cm 




nidi = 2.75 


[104 K] 


0.9 


Wr 


( 2 ) 


where T is the temperature of the ionized gas, iie is the lo 
cal electron density, and L is the total p ath length. Assuming 
T = 104 K at high Galactic latitudes (e.g. 


Sun et al. 


2008 


Aka- 


jhori et ari|2013[ and references therein), we derive the related 
emission measure of EM - \.l ± 0.3 pc cm“® across the 3C196 
field. This corresponds to a large-scale mean electron density of 
(mc) ^ 0.03 cm“^ over an assumed path length of 1.5 kpc. Its 
spatial variation across the 3C196 field is <j(n^) ^ 0.003 cm“^. 
In our calculation, we do not take observed vertical temperature 
gradient in the diffuse ionized gas into account ( [Reynolds et al. 
[1999[|Haffner et al.[2009| l, as the gradient is small and our results 
are not sensitive to it. 

We can now use the equation for rotation measure. 


(RM) 
[rad m^^] 


= 0.81 


I 


<»e) (^ll) 

0 [cm-3] [pG] 


dl. 


(3) 


and estimate cr. 


■(Bii) according to 


cr 




^■(RM) \ I {RM}o-(n,) 


0.81 {n,}L) \o.8l {n,fL 


(4) 


Schnitzeler 


(20101 


Rc(®) = 




(5) 


where sums are taken over all pixels of an image S jj at Eara- 
day depth <6, where Sij >1.5 mJy PSE^' RMSE~^ (at least two 
times above the noise). The vector rij contains the coordinates 
of the pixels. 

We find that the centroid of emission changes systematically 
by ~ 0.8° per rad m“^ from west to east and by 0.3° per rad 


1 2003 [ l survey, which means that given Ho- emission is at 1° resolution. 
The WHAM survey provides velocity resolved spectra, but for this work 
we use only integrated emission. 


from north to south. Erom Earaday depth H-1 to -1-2.75 rad 
the polarized structures have moved by 1.4° towards the east and 
by 0.6° towards the south. This is indeed an indication of a large- 
scale gradient in RM. 

Based on our previous discussion that the ct^rm) term in Eq.|^ 
dominates over the cr^n,) term, the large-scale gradient that we 
observe seems to be a gradient of across the 3C196 field. 
However, the distribution of the thermal plasma is known to be 
extremely inhomogeneous with a very small filling factor for 
the Warm Ionized Medium (WIM, e.g. Berkhuijsen et al.|2006 


In our case, the first term in Eq. dominates over the second 
term. Using L = 1.5 kpc and (RM) = 1 rad m“^, we get varia¬ 
tions of the mean line-of-sight magnetic field component across 
the 3C196 field to be - 0.1 pG. This is in agreement with 
an upper limit, ^ 0-4 juG, derived by ^ “ 

using the rotation measure of NVSS sources.' 

If B||, rig, or both show a large-scale gradient across the sky, 
this should be seen in observed RM cubes. Polarized structures 
appear to move in the direction of the gradient as we go through 
different Earaday depths. We see this behaviour in our RM cubes 
from H-1 to h-2.75 radm^^. Polarized emission moves mostly 
from the west to the east (see Pig. 0&|6). To quantify this, we 
calculate the centroid of emission at each Earaday depth as 


[Gaensler et al.|2008), which is mostly responsible for the Para 
day rotation (e.g. Heiles & Haverkorn 2012| l. The estimated tig 
from the [Pinkbeiner ( [2003[ l Ha map is just an average of the true 
distribution along a hne-of-sight. Therefore, it is possible that 
inhomogeneity of rig also contributes to the observed gradient in 
Earaday depth. 

Pinally, By across the 3C196 field varies both in strength 
and in orientation. This is seen through the presence of ob¬ 
served structures both at positive and negative Earaday depths 
(see Pig.|^. Positive RM values are associated with By orientated 
towards us and negative RM values with By orientated away form 
us. This does not necessarily mean that By has the same orien¬ 
tation along the line of sight. The field can reverse, sometimes 
even multiple times. 

Direct evidence for field reversal(s) can be seen in the Para- 
day spectrum, where structures are found at both positive and 
negative depths. In our RM cubes there are ~ 5% of Earaday 
spectra that we directly associate with field reversal(s) (for an 
example see Pig.|^. However, spectra with structures, which are 
only on one side of the Earaday spectrum, do not exclude field 
reversals. If two regions with opposite magnetic field orienta¬ 
tions change the polarization angle of incoming emission, but 
they do not themselves emit, and there is no additional emission 
between the two, the polarization angle of incoming emission is 
first rotated in one direction and then derotated in the other. The 
observed polarization angle is then a net effect of Earaday rota¬ 
tions along these two regions (for an illustration, see Pig. 2 in 
Brentjens & de Bruyn|2005). In the Earaday spectrum, this only 


manifests as one structure, despite the fact that in reality we have 
two physically separated regions. 

5.2. Physical picture and relative distribution of intervening 
magneto-ionic medium 

The synchrotron emission of our Galaxy is known to have many 
components: the disk, a thick disk, and possibly a halo (e.g. 


Subrahmanyan & Cowsik 2013 and references therein). At high 


Galactic latitudes the bulk of the emission comes from the thick 
disk with a width of about 1-1.5 kpc. The polarized emission 
builds up over the same path length, but along the way it gets de¬ 
polarized by the magneto-ionic ISM. Depolarization is stronger 
at lower frequencies. A typical observed polarization at LOPAR 
frequencies is about 1 - 4% ( [lacobelli et al. 2013[ Jelic et al. 
[2014[ l, while intrinsically synchrotron radiation is polarized at 
the 70% level. Hence only a small fraction of the synchrotron 
emission of our Galaxy is observed in polarization. 

A possible geometric configuration and distribution of the 
random and uniform component of the magnetic field and 
magneto-ionic medium are very difficult to get directly from ob¬ 
servations. Analysis of observations in Earaday depth helps to 
constrain the relative distribution of structures. However, the dis¬ 
tribution in Earaday depth does not necessarily reflect the true 


physical distribution (see discussion in Brentjens & de Bruyn 
[2005[l. Magnetic field reversals and inhomogeneous distribution 
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of regions (i) that only emit; (ii) with only Faraday rotation; and 
(iii) with a combination of both, give a rise to a discrepancy be¬ 
tween the two distributions. Full radiative transfer simulations 
are needed to disentangle this problem. 

To visualize the complexity of the problem, we have con¬ 
structed a simple physical model (see Fig. [^, which we believe 
contains the essential ingredients to explore the conditions in 
the ISM that are required to explain observed structures in the 
3C196 field. Observed features are organized in four distinc¬ 
tive groups; (i) filaments A and B: relatively straight filamentary 
structures running from north to south at Faraday depths around 
-1-0.5 rad m“^ and parallel to the Galactic plane; (ii) diffuse emis¬ 
sion BGl; prominent diffuse background emission at Faraday 
depths from -i-l.O to h- 4.5 rad m“^, showing a large-scale gra¬ 
dient; (iii) diffuse emission BG2: very faint patchy emission at 
Faraday depths from -1-4.5 to h- 7.5 rad m“^ and at negative Fara¬ 
day depths from -3 to -0.5 rad m“^ (iv) a horizontal structure 
located in the upper north of the field; and (v) a “triangular” fea¬ 
ture C in the south-west. 


We now describe and discuss our crude model. At first, we 
assume that dominant orientation of By is towards us. This is 
consistent with the observed, mostly positive, values of Fara¬ 
day depth. Based on Eq. assu ming a path length of 1.5 kpc, 
taking estimated values from Sec. |5.l| of {n,,) ^ 0.03 cm~^ and 
{B|l) ^ 0.3 and including a vertical gradient in electron den¬ 
sity 1^ , we expect to see polarized emission at Faraday depths 
up to d) < -1-10 rad m“^. This is in an agreement with our ob¬ 
servations, where we have detected faint patchy structures up 
to -1-8 rad m“^. Simulated RM profiles towards the 3C196 field 
(/ = -1-171, b - h- 33°; see Fig. 12 in [Sun et al. 2008| l, based 
on the regular double-torus halo field and disk field models de¬ 
scribed in Sun et alT] ( 2008| l, also agree with our observations. 
These profiles predict RM values up to < H-10 rad m“^. 

At the negative side of the Faraday spectrum, structures are 
only observed up to -3 rad m“^. Under the same assumptions, 
but with the magnetic field pointing away from us, the thick¬ 
ness of the negative Faraday rotating region appears to be much 
smaller than the positive: Eq. |^gives L ^ 250 pc. 


A relative distribution of the observed features is as follows. 
Eilaments seem to be displaced from the background diffuse 
emission BGl by only ~ 1.5 rad m“^. Their brightness is com¬ 
parable to the background emission that shows a very uniform 
large-scale structure of polarization angle on both sides of the 
filaments (see high-resolution images in Fig. [^. We therefore 
conclude that filaments are non-emitting magnetized plasma lo¬ 
cated in front of the prominent background BGl. 

The Faraday depth of the filaments is smaller than that of the 
background emission. This implies By with an orientation in the 
opposite direction to the background. Parts of the background 
emission are derotated within the filaments, giving a deficit in 
Faraday depth compared to the other parts of the background 
emission. 


Prominent background emission BGl is closer to us than 
the fainter background emission BG2, and is probably a very 
nearby component. BGl emission might be associate with the 
synchrotron emission coming somewhere from the edge of the 
Local Bubble, which is estimated to have a radius of approxi¬ 
mately 200 pc (see e.g. |Sun et aL]|2008|l. To be consistent with 


’ The vertical electron density gradient can be described as n^iz) = 
where is a mid-plane electron density and h is a scale height 
(e.g. |Haffner et ^.|2009[|Schnitzeler|2012[ and references therein). In 
our calculation, we normalize it to give {n^) ^ 0.03 cm“^. 



Fig. 8.A cartoon of the possible layout of polarized emission and mag¬ 
netic fields towards the 3C196 field. Solid arrows indicate orientation 
and strength of By. For display purposes, the filamentary structures are 
shown along the viewing direction in this sideways projection. 


our observations, in front of this emission there must be a non¬ 
emitting magnetized plasma, with By orientated towards us. 

The same orientation of By we also expect for a region asso¬ 
ciated with the fainter emission BG2. This emission is the most 
distant component and probably reflects the global emission of 
the thick disk. Along the way, two types of depolarization occur 
(for an overview, [Haverkorn et al.|2004| and references therein), 
resulting in the observed very faint and patchy emission. The 
first is wavelength-independent depolarization, due to turbulent 
magnetic fields in the ISM, and the second is due to differential 
Faraday rotation occurring in regions where emitting and Fara¬ 
day rotating areas are mixed. 

The horizontal structure in the north, at negative Faraday 
depths, does not have any positive Faraday depth components 
along the line of sight and, therefore, it is a region with a dom¬ 
inant By component orientated away from us. This horizontal 
structure also seems to be spatially a part of the large-scale gra¬ 
dient observed in BGl. If this is true, in front of the synchrotron 
emitting region BGl, there is one layer of a non-emitting magne¬ 
tized plasma with By orientated towards us and one layer with By 
orientated away from us. The second region should have a large- 
scale gradient in By. Since the strength of By along the gradient is 
different, some part of the background emission appears at posi¬ 
tive and some at negative Faraday depths. Then, in this case, the 
observed filamentary structures and the horizontal structure in 
the north reflect the regions with an excess in thermal electron 
density and/or magnetic field strength. However, the filaments 
should be located in the front of all other components since they 
cast shadows both across emission observed at positive and neg¬ 
ative Faraday depths. 

Finally, Faraday spectrum at a location of the negative tri¬ 
angular structure C peaks both at positive and negative Fara¬ 
day depths. Positive Faraday components are associated with 
the background emission BGl and BG2. The negative triangular 
component is probably associate with an additional emitting re¬ 
gion, located somewhere in front of the background emission 
BGl and BG2. We tried to incorporate all of these into our 
model. A cartoon of this model is shown in Fig. 

5.3. Possible associations with the filamentary structure A 
and its proper motion 

The elongated straight filament A running in the north-to-south 
direction is displaced by ~ 1.5 rad m“^ in Faraday depth relative 
to the surrounding background emission, suggesting a magneto- 
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Fig. 9. Examples of higher resolution images of the 3C196 field in polarized intensity (PI) and polarization angle (PHI) given at Faraday depths 
of +1.0, +1.5, and +2.0 rad Images show the central 2.5° x 2.5° of the field. During the imaging we used baselines between 10 and 4000 
wavelengths, providing a frequency-independent resolution of about 1 arcmin. 


ionic structure located in front of the bulk of the emission. It is 
interesting to compare it to other linear ISM features that have 
been reported in the literature. iMcCullough & Benjamin (2001 1 
discovered an extremely narrow, linear, ionized filament in the 
ISM, which resembles some characteristics of our filament at 
Faraday depth +0.5 rad m“^. Their filament is 2.5° long, 20" 
wide and has an Ho- surface brightness of Ih„ ~ 0.5 R. It is 
located at a comparable Galactic latitude {b — +38.5°) to the 
3C196 field. IMcCullough & Benjaminj (|2001|l argued that the 
most probable origin of the filament is an ionized trail left by 
photoionization from a star or compact object. As an alternative, 
but unlikely scenario, they mentioned an extremely low density 
nearby stellar jet, or an unusually linear filament associated with 
a large-scale nearby bubble, or an ionized trail left by mechanical 
input from a star or compact object. They also argued that this 
kind of filament would show a displacement in Faraday depth 
against the background radio emission, regardless of its true na¬ 
ture of origin. 

Unfortunately, up to now there are no observations in the ra¬ 
dio that can support this prediction. This is similarly true for our 
filament. Available data in Ha of 3C196 field do not have the 
angular resolution needed to resolve the filament. For the pur¬ 
pose of further discussion, we assume that our filament has the 
Ha surface brightness of the [McCullough & Benjamin pOOlj l 
filament, and then estimate the thermal electron density of the 
filament and its distance. 

Combining Eq.|^&[^ we find that the thermal electron den¬ 
sity of the filament is 


rip EM 

= O.8IB1 


[rad m 


-2n 


[cm“2] 


RM [/zG][pc cm-6] ’ 


(6) 
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where to calculate rip we need to assume the strength of B\\. Once 
that rip is known, one can get dl from Eq.|^or[^and then estimate 
a distance, d, to the filament based on its angular diameter, Dg - 
9.4', and assuming axial symmetry. 

Taking Ih„ - 0.5 R and RM =1.5 rad m“^ with the results 
from the previous sections, T — lO'^ K and - 0.3 fiG, we 
get EM - 1.37 pc cm“®, rip - 0.2 cm“^, dl - 28 pc, and d - 
20 kpc. The calculated distance is obviously not realistic. The 
filament cannot be located behind the background emission (see 
discussion in Sec. |5.2| ). This means that some of our assumptions 
in the calculation are not correct. 

Eirst, it is possible that our filamentary structure is not a fil¬ 
ament but rather a sheet-like structure observed edge-on. Axial 
symmetry in that case is not valid, and the distance based on 
this assumption would not be valid either. Second, if we keep 
the assumption about the filamentary nature, then either B\\ or tip 
or both must be larger. If we assume that our filament is located 
somewhere within the Local Bubble, i.e. d < 200 pc, this implies 
that dl < 0.3 pc and npB\\ > 6.2 cmT^jiG. A filament is a region 
showing an excess in rip and/or B\\ compared to its surroundings. 

Eor a proper understanding of the nature of the filament, it is 
crucial to know its distance and proper motion. We assume that 
the filament is at 50 pc distance and has a transverse velocity 
of 50 km s *. It should then move 50 ■ 10“® pc year ’, which is 
equivalent to 0.2" year '. This causes a parallax with an ampli¬ 
tude of about 0.03" over the course of one year. Obviously it is 
easier to measure the filment’s proper motion than its distance. 
We therefore try to measure the first using our 1 arcmin reso¬ 
lution images (see Eig. |^, observed one year apart (L80508 & 
L192832 observations). 
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The edges of the filament are well defined by the depolariza¬ 
tion canals seen so clearly in the upper left image in Fig. 0 We 
use these canals to create a mask and isolate the filament from 
its surrounding background emission. We repeat this for the two 
observations and then determine if the position of the filament 
has changed on the timescale of one year. We find no change in 
position between the L80508 & L192832 observations. The es¬ 
timated error on our current measurement is ~ 1". Assuming a 
typical transverse velocity of 50 km s“', we can set a lower limit 
to the distance of > lOpc. 

To achieve the precision needed to measure a real proper mo¬ 
tion, we can reduce the noise by averaging the RM cubes from 
many nights in each of the two seasons. However, for this to 
work for a very small fraction of the width of these filaments, it 
is essential to compare the structures at exactly the same Faraday 
depth. The RM cubes must therefore be aligned in RM to much 
better accuracy than provided by the CODE ionospheric data. 
Recently, Brentjens et al. (in prep.) developed a new method, 
which yields the required accuracy limited only by the available 
signal within the field of view. The method cross-correlates the 
Galactic polarized emission in Faraday depth as a function of 
time, both within each night and between nights and season. Ini¬ 
tial results suggest that the ionospheric Faraday rotation can be 
aligned to levels on the order of 0.01 rad m~^, i.e. a very small 
fraction of the width of the RMSF The results of this analysis 
will be presented in a future paper. 


5.4. The system of depolarization canals 

Figure|^shows many organized patterns of linear structures with 
no or very weak polarized signal. Many of these resemble the 


depolarization canals first discussed and studied by Haverkorn 
et al. (|2000| |2004 1 . We also refer to these structures as canals. 
Haverkorn et al.| ( 2004| l described two characteristics of these 
depolarization canals. They are one PSF wide and the polariza¬ 
tion angle changes across the canal by 90°. Mechanisms that cre¬ 
ate them are as follows. The canals are at the boundary between 
two regions, which each have similar polarized intensities but a 
difference in polarization angle of A0 = (n H- 1/2)180° (where 
n = 0,1,2,...). This causes almost total depolarization through 
vector addition, and the resulting canal is by definition one PSF 
wide. The other mechanism is connected to differential Faraday 
rotation. A linear gradient in RM produces a long, narrow depo¬ 
larization canal at a certain value RMc, where 2RMc4? - nn. 
Across every null in the sine-function, the polarization angle 
changes by 90°. A third possibility is beam depolarization due 
to the presence of turbulent magneto-ionic medium, which is ex¬ 
tremely small-scale, in front of a polarized background. To ex¬ 
plore this possibility, we also made a RM cube with a spatial 
resolution of 12". At this scale the canals appear resolved, but 
they do not show an increased polarization signal. We therefore 
consider this possibility unlikely as a cause for the creation of 
canals. 

We believe that most of the canals are probably of the first 
kind, that is they are due to beam depolarization due to disconti¬ 
nuity in the polarization angle. The best example are the canals 
associated with the edges of the prominent filament discussed 
previously and the spaghetti-like structure located to the east of 
it (see Fig.[^. They are indeed one PSF wide, the intensity of the 
emission on either side of the canal is comparable, and there is 
a 90° polarization angle change when crossing the canal. Adopt¬ 
ing this mechanism for the creation of a canal there is, however, 
a clear prediction: the polarized emission cannot go to zero at 
all frequencies that were used in the creation of the RM cube. 


which by necessity spans a wide range of frequencies (in our 
case from 115-175 MHz). The Stokes Q and U signals at the 
edges of this frequency range should therefore increase again 
within the canals. Unfortunately the signal-to-noise ratio in our 
data at this stage is not yet sufficient to confirm this prediction. 
We stress that the residual polarized signal within the canals also 
carries valuable information about the polarization signal emit¬ 
ted between the location of the filament and the observer at any 
Faraday depth. In this picture the deepest canals must be clos¬ 
est to the observer. We will return to this interesting aspect in a 
future paper when we have aligned and analysed, the RM cubes 
from the many nights we have accumulated on this field. 

There are also some canals in our observations that are 
broader than one PSF. One of these is located in the south-west 
part of the image, along the boundary between a triangular re¬ 
gion showing emission at negative Faraday depth and surround¬ 
ing emission at positive depths (see Fig. |^. This type of canal 
must be either a region in the ISM showing no emission, all 
along the line of sight, or a region with an extremely turbulent 
magneto-ionic medium, causing wavelength-independent beam 
depolarization. 


5.5. Straightness of observed structures and canals 


Without regard to the origin of observed structures and canals, 
a real puzzle is their extraordinary straightness. Their axial ratio 
is larger than 100:1. It is possible that they are the results of a 
fortunate projection of a three-dimensional morphology of the 
magnetic field (e.g. the folding of magnetic sheets and loops) 
and the ISM. 

Recent magneto-hydrodynamical simulations of the ISM, by 
Choi & Stone ( 2012| l, showed that thermal conduction plays an 
important role in shaping the geometry of structures formed by 
the thermal instabilities by stabilizing small scales and limiting 
the size of the smallest condensates. If the magnetic field is uni¬ 
form and weak (~ 3 nG), the heat is conducted anisotropically, 
primary along magnetic field lines. As a result of this, the ther¬ 
mal instabilities saturate into very long thin filaments of dense 
gas aligned with the magnetic field. In the linear regime, this can 
be explained by the variation of the Field length ( |Field|[T965] l 
with respect to the direction of magnetic field. Isobaric perturba¬ 
tions with wavelength smaller than the Field length do not grow 
at all. Perpendicular to the magnetic field, the Field length is very 
small, so very short wavelength perturbations are unstable. Par¬ 
allel to the magnetic field, the Field length is much longer, so 
only longer wavelength modes grow. As the perturbations grow 
non-linear, conduction tends to enforce isothermality along mag¬ 
netic field lines, leading to long filaments. If the magnetic field is 
uniform and strong (~ 3 pG), only motions along field lines are 
allowed, and, moreover, magnetic pressure provides some sup¬ 
port in dense regions, resulting in fragmented filaments aligned 
with the field lines (for more details we refer to |Choi & Stone| 


2012| l. Therefore, it is possible that filamentary structures ob¬ 
served in our data were shaped through thermal instabilities. To 
get better statistics on observed structures/filaments, and to study 
their distribution and orientation with regard to the regular mag¬ 
netic field component, we plan to survey a larger area of the sky. 
This will then allow us to directly associate observed structures 
with the results of magneto-hydrodynamical simulations of the 
ISM. 

Another explanation of observed linear structures can be 
connected to interactions of close-by moving stars with the ISM, 
as discussed by McCullough & Benjamin ( 2001| l. Bow shocks 
around stars and neutron stars (pulsars) can leave long and nar- 
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row trails in the ISM, because of their large speeds. Or, trails can 
be left by an ionizing source of sufficiently low ionizing lumi¬ 
nosity that its Stromgren radius is small compared with the dis¬ 
tance it travels in a recombination time. However, to make such 
trails long in angular size, of a few degrees, these sources should 
be rather close, less than about 100 pc, and located somewhere 
within the Local Bubble. We will leave a more detailed study of 
this possibility for future work. 

6. WSRT observations from 312-381 MHz 

The 3C196 field was also observed with the WSRT at 350 MHz 
in 2013 and 2014. These observations form part of a large-area 
mosaic survey to complement LOFAR polarization studies. Thus 
far we have finished the reduction and analysis of the part, which 
overlaps with the LOFAR field of view so we can make a com¬ 
parison with the features seen at LOFAR frequencies. 

Figure [T^ shows the WSRT 350 MHz image in polarized in¬ 
tensity at Faraday depth 0 rad m“^. The RMSF in the WSRT 350 
MHz data (dO =11 rad m“^) is an order of magnitude less than 
in the LOFAR data. Hence, almost all Faraday depths that show 
emission in the LOFAR data are captured in this one frame in the 
WSRT data. The linear filament A passing over 3C196 is clearly 
detected. It has a surface brightness of ~ 7 mJy PSF ' RMSF^^ 
This is the first time that we see some common morphology in 
the RM cubes made in the LOFAR 115 - 175 MHz and WSRT 
315-380 MHz bands. Diffuse emission surrounding the filament 
is not as prominent (< 3 mJy PSF'' RMSF'*) as observed with 
LOFAR (see Fig. Q and its morphology is very patchy. 

In addition to the filament A passing over 3C196, we de¬ 
tected yet another straight filamentary structure (marked with 
D in Fig. m at the WSRT frequencies. It is located in 
the south-east part of the image. Its surface brightness is ~ 
8 mJy PSF'* RMSF * and is probably associated with emission 
at -1-2.5 rad m'^ in the LOFAR RM cubes. The two filaments in 
the WSRT data seem to have similar orientation. 

Bright horizontal emission (~ 15 mJy PSF * RMSF'*) lo¬ 
cated to the north of the filaments in the WSRT data coincides 
with the region showing weak emission and horizontal depo¬ 
larization canals at LOFAR frequencies, just below horizontal 
emission in the north at negative Faraday depths. Some parts of 
this bright emission in the WSRT data are probably Faraday- 
thick structures at LOFAR frequencies and not detectable with 
LOFAR (Affi > ^ffiLOFAR)- 

The different observed morphology of the diffuse emission 
surrounding the filament in the two different frequency bands 
may be attributed to the poor resolution in Faraday depth at 350 
MHz. To test this, we generated an RM cube at LOFAR fre¬ 
quencies using an RMSF that has a resolution of (5®35o mhz ~ 
10 iSOiso MHz- It is possible that multiple Faraday thin structures 
(AO < i5035omhz) detected in the LOFAR low-frequency images 
decorrelate when we observe them with a much broader RMSF. 
We find that ~ 50 % of the flux is depolarized, but the overall 
morphology of emission is not much different from the origi¬ 
nal LOFAR RM cubes. Thus, the lack of prominent background 
emission detected with LOFAR cannot be explained by the much 
broader RMSF at 350 MHz. 

Another possible explanation for the lack of emission in the 
WSRT 350 MHz data may be related to the missing large-scale 
structures. The shortest baseline of the WSRT is 36 m, so at 350 
MHz structure on angular scales larger than about a degree is 
not adequately measured. As we discussed in Sec. the syn¬ 
chrotron radiation is emitted on large scales, but its polarization 
is altered on smaller scales by Faraday effects. Faraday rotation 


WSRT 350 MHz, PI: <i> = 0 rad/m^ 



128 126 124 122 120 

a (J2000) 

Fig. 10. An image of the 3C196 field in polarized intensity (PI) observed 
with the 28 WSRT mosaic pointings at 350 MHz (based on preliminary 
calibration and analysis). In making this image, we used only the cen¬ 
tral 9 of the 28 observed pointings. The image is given at Faraday depth 
0 rad m'^ and has not been corrected for polarization noise bias. It mea¬ 
sures 6.9° X 6.9° in size and the PSF is 2.8 x 3.8'. 


is much stronger at lower radio frequencies. We thus expect the 
missing large-scale emission in polarization to be more promi¬ 
nent at higher frequencies. This was discussed first by |Haverkorn| 
et al.| ( 2004| . For the case of a uniformly polarized background 
propagating through a small-scale Faraday screen, the expected 
offset for Stokes Q and U depends on the wavelength. A, and 
the width of the Faraday depth distribution, cr^, as e~ i . In 
WSRT 350 MHz observations, this condition is not well satis¬ 
fied and one needs to take it into account during interpretation of 
the data (e.g. Haverkorn et al.|2004 [Schnitzeler et al.|2009) l. 

To correct for these undetectable large-scale structures, one 
needs first to observe the same region at the same frequen¬ 
cies with a single-dish telescope with absolute intensity scaling. 
Then, these large-scale data can be added to the interferometer 
data ( Uyaniker et al.|199^|Stanimirovic|2002| l. However, for the 
WSRT 350 MHz observations this is not possible. We can only 
estimate it based on other existing absolutely calibrated obser¬ 
vations at higher frequencies, e.g. |Berkhuijsen & Brouw ( |1963| l 
map of polarized emission at 408 MHz. This frequency is close 
enough to 350 MHz to allow a comparison, although the polar¬ 
ized intensity at 408 MHz is expected to be slightly higher due 
to the larger polarization horizon. 

The polarized brightness temperature at 408 MHz at the lo¬ 
cation of the 3C196 field is 2.1 + 0.4 K ( [Berkhuijsen & Brouw| 
1963| l. Using a power-law spectral index of -2.5, this corre¬ 
sponds to ~ 3 K at 350 MHz. If we now smooth our 350 MHz 
data to the 2° FWHM of Berkhuijsen & Brouw| ( |1963| l we get 
0.08 K. This means that indeed we have a missing large-scale 
component in polarized intensity. 


7. Summary and conclusions 

We have presented results from LOFAR HBA observations of a 
field centred on the 3C196 source, taken as part of the LOFAR- 
EoR project. We applied RM synthesis to the data to reveal a 
very rich morphology of polarized emission. Detected structures 
are spread over a wide range of Faraday depths, ranging from -3 


Article number, page 14 of 16 






































V. Jelic et al.: Linear polarization structures in the 3C196 field 


to +8 rad tn“^. Their average brightness temperature is 5 - 15 K, 
which is comparable to the emission in the ELAIS-Nl field and 
also to a field observed with LOFAR at high Galactic latitudes 
Pelic et al.|20T^ . 

The most interesting features detected in polarization are 
(i) a straight and long filamentary structure at Faraday depth 
+0.5 rad nT^, parallel to the Galactic plane and (ii) a system 
of linear depolarization canals. The first is probably an ionized 
structure in the ISM located very close by, somewhere within 
the Local Bubble. It is displaced by 1.5 rad m“^ from the sur¬ 
rounding background emission. This implies that, within the fil¬ 
ament, the thermal electron density and magnetic field must sat¬ 
isfy neB\\ > 6.2 cmr^jiG. With the current accuracy of the posi¬ 
tion measurement of 1", we can set a lower limit to its distance 
of 10 pc, assuming a transverse velocity of 50 kms^'. 

The same filament was also observed with the WSRT at 
350 MHz. This is the first time that we see some common 
morphology in the RM cubes made in two different frequency 
regimes. Interestingly, in the south-east part of the WSRT mo¬ 
saic there is another filament. This filament is probably associ¬ 
ated with emission at +2.5 rad m“^ in the LOFAR RM cubes. Its 
orientation is the same as the filament located in the centre of 
the image. This is a tentative indication that these two filaments 
might be a part of the same large-scale structure in the ISM. 

A follow up study of the filaments discussed in this paper 
showed their strong correlation with the orientation of the sky 
projected magnetic field component as probed by the Planck 
maps of the dust emission in polarization ( |Zaroubi et al.|2015| l. 
Dust particles are usually embedded within neutral hydrogen 
(Hi) gas, as seen through Hi/dust correlation at high Galactic lati¬ 
tudes ( |Planck Collaboration et al.|2014) l. Therefore, high angular 
resolution Hi data of 3C196 field can be very usefully in provid¬ 
ing potential kinematic information about the features studied in 
this paper. This is also supported by [Clark et'af] ( |2014[ ) anal¬ 
ysis of Hi data from the Galactic Arecibo L-Band Feed Array 
Hi (GALFA-Hi) Survey and the Parkes Galactic All Sky Sur¬ 
vey (GASS), which revealed many slender, linear Hi features 
(“fibers”) that extent for many degrees at high Galactic latitude. 
These fibers are orientated along the interstellar magnetic field as 
probed by starlight polarization and are most likely a component 
of the Local Bubble wall. 

Most of the depolarization canals detected in our LOFAR 
data are due to beam depolarization. However, some of them 
might be associated with regions in the ISM that show a very tur¬ 
bulent magneto-ionic medium, causing wavelength-independent 
beam depolarization. Without regard to their origin, the straight¬ 
ness of these canals is a real puzzle. Are they a projection effect 
of complicated morphology in the magnetic field and the ISM, 
or are they associated with fast-moving close-by stars interacting 
with the ISM? This needs to be investigated further. Planned ob¬ 
servations covering a much larger area of this region, along with 
observations at other wavelengths, will help towards this goal. 

Results presented in this paper are also of concern for epoch 
of reionization experiments. The instrumentally polarized re¬ 
sponse of an interferometer needs to be calibrated to a small frac¬ 
tion of a percent to limit leakage of polarization signals to levels 
of a few mK. Otherwise, the leakage can contaminate the cos¬ 
mological 21-cm signal in total intensity (e.g. [Jelic et al.|20T0j l. 
The contamination depends strongly on the brightness temper¬ 
ature of polarized emission, its morphology, and its spread in 
Faraday depth. In the case of the polarized emission detected 
in the 3C196 field, which is limited to sm all Faraday depths 


(|<h| < 10 rad m ^), this is not a big concern. Asad et al. ( j2015 i 
showed that in the spherical power spectrum the leakage power 


of emission presented in this paper is lower than the cosmolog¬ 
ical signal at k < 0.3 Mpc Moreover, there is a window in 
the cylindrical power spectrum where the cosmological signal 
dominates over the leakage. 
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